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Exponential  suppression  of  thermal  conductance  using  coherent  transport  and  heterostructures 
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We  consider  coherent  thermal  conductance  through  multilayer  photonic  crystal  heterostructures,  consisting 
of  a  series  of  cascaded  nonidentical  photonic  crystals.  We  show  that  thermal  conductance  can  be  suppressed 
exponentially  with  the  number  of  cascaded  crystals  due  to  the  mismatch  between  photonic  hands  of  ah  crystals 
in  the  heterostructure. 
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Coherent  thermal  transport,  using  thermal  channels  with 
length  smaller  than  the  mean-free  path  of  the  thermal  carri¬ 
ers,  is  important  for  fundamental  study  of  thermal  processes 
and  for  new  device  opportunities  in  thermal  management.'  " 
In  this  Brief  Report,  we  consider  a  coherent  thermal  channel 
formed  by  cascading  a  series  of  nonidentical  photonic  crys¬ 
tals.  We  show  that  its  thermal  conductance  can  be  suppressed 
exponentially  with  respect  to  the  channel  length  1.  This  is 
fundamentally  different  from  the  incoherent  process,  where 
the  conductance  typically  decreases  as  1//.'^  We  therefore 
demonstrate  that  coherent  processes  can  be  very  effective  in 
suppressing  thermal  conductance.  This  result  also  indicates 
that  an  aperiodic  coherent  thermal  channel  is  qualitatively 
different  from  all  previously  considered  coherent  thermal 
channels'  *"  including  periodic  multilayer  photonic  crystals 
as  we  previously  considered,"’*"  where  the  intrinsic  thermal 
conductance  of  the  channels  were  all  independent  of  the 
channel  length  /. 

As  a  concrete  implementation,  we  utilize  the  concept  of 
photonic  crystal  heterostructure, *^  *^  and  consider  multilayer 
systems  consisting  of  a  total  of  N  different  photonic  crystals, 
as  illustrated  in  Fig.  1(a).  All  crystals  consist  of  alternate 
periodic  layers  of  vacuum  and  dielectric.  The  use  of  vacuum 
ensures  that  heat  transfer  is  carried  only  by  photons.  The 
dielectric  can  be  silicon  (n  =  nsi=3.42).  For  photon  frequen¬ 
cies  within  the  blackbody  spectrum  at  room  temperature,  sili¬ 
con  has  very  little  dispersion  and  dissipation,  with  typical 
attenuation  length  exceeding  millimeter.'"  Consequently 
photonic  thermal  transport  should  be  coherent  across  the  en¬ 
tire  structure."’*" 

For  the  structure  shown  in  Fig.  1(a),  we  assume  that  all 
crystals  have  the  same  period  a.  The  mth  crystal  has  a  di¬ 
electric  layer  thickness  d^,n  and  a  vacuum  layer  thickness 
At  either  ends  of  the  structure,  we  have  two  semi¬ 
infinite  photonic  crystals.  In  between  these  two  ends,  there 
can  be  a  series  of  crystals  cascaded  together,  each  having  the 
same  number  of  periods  N-p  and  hence  the  same  total  thick¬ 
ness.  (Thus  the  length  of  the  channel  is  proportional  to 
N-2.)  We  also  assume  that  the  two  crystals  at  the  ends  are 
maintained  at  temperatures  of  T  and  T+dT,  respectively.  The 
three-dimensional  (3D)  thermal  conductance  per  unit  area  is 
then  defined  as  G^p,{T)  =  dQ{T)/dT,  where  dQ{T)  is  the  heat 
flux  per  unit  area. 

In  such  multilayer  structures,  each  photon  state  is  charac¬ 
terized  by  three  parameters:  frequency  w,  wave  number  k\\  in 
the  direction  parallel  to  the  layers,  and  polarization  a=s,p. 
The  s  and  p  polarizations  have  electric  and  magnetic  fields 


parallel  to  the  layers,  respectively.*"  Summing  over  all  pho¬ 
ton  states,  we  have*" 


G3d(T)  =  E  f  ^G(7’,k,|,a’).  (1) 
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Here 


G(r,k||,(r)  =  kB 
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©(o),  k||,  O’) 


(2) 


is  the  thermal  conductance  of  a  one-dimensional  thermal 
channel  of  a  given  parallel  wave  number  ky  and  polarization 
cr.  kp  is  the  Boltzmann  constant  and  h  =  hl{2TT)  is  the  re¬ 
duced  Planck  constant.  The  factor  0(a),k||,cr)  measures  the 
contribution  from  photon  states  at  {co,k\\,a-).  For  a  single 
crystal  (i.e.,  N=  1),  @(&i,k||,  cr)=  1  if  the  frequency  co  lies  in  a 
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FIG.  I.  (a)  A  photonic  crystal  heterostructure  with  four  crystals. 
Each  crystal  has  alternate  silicon  and  vacuum  layers,  and  a  period  a. 
The  two  crystals  at  either  ends  are  assumed  to  he  infinite,  (b)  In  the 
top  four  panels,  gray  regions  correspond  to  frequency  ranges  where 
the  corresponding  crystal,  assumed  to  be  infinite,  has  a  propagating 
state  along  the  normal  incident  direction.  White  regions  are  the 
hand  gaps,  is  the  thickness  of  the  silicon  layers.  In  the  bottom 
panel,  the  hlack  regions  represent  frequency  ranges  where  there  are 
photonic  bands  in  all  crystals. 
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FIG.  2.  (Color  online)  (a)  Normalized  one-dimensional  thermal 
conductance  G{T)/Gq(T)  as  a  function  of  temperature  T  for  hetero¬ 
structures  with  different  number  of  crystals  N.  For  each  N,  different 
curves  represent  structures  of  different  thicknesses  of  silicon  and 
vacuum  d^.  The  brown  (dark  gray)  dashed  curve  shows  a  commen¬ 
surate  case  of  N=2,  where  iisi=0.7851a  and  rfs2  =  0-2262a  such  that 
”si‘^s2/(®“‘^s2)=  1-  (b)  High-T  conductance  for  heterostructures  of 
different  number  of  periods  Np  at  the  intermediate  crystals.  Each 
line  corresponds  to  structures  with  the  same  number  of  crystals  N. 
3£A£8.  Data  points  exist  only  at  integer  Np.  (c)  High-T  conduc¬ 
tance  for  a  total  of  15  000  heterostructures  with  randomly  generated 
thicknesses  of  silicon  layers.  0.05a£d^£^0.95a  for  all  crystals. 

photonic  band,  and  @(w,A:||,  cr)  =  0  if  the  frequency  w  lies  in  a 
band  gap  region. For  the  heterostructure  [Fig.  1(a)],  where 
there  are  at  least  two  crystals  (i.e.,  A&2),  0(w,^||,cr)  van¬ 
ishes  unless,  at  a  given  (w,A:||,cr),  propagating  eigenmodes 
exist  for  the  semi-infinite  crystals  at  both  ends.  In  such  a 
case,  @(<u,A:||, cr)  is  the  power  transmission  coefficient 
through  the  entire  structure  when  an  incident  wave  is  one  of 
these  eigenmodes,  and  can  be  determined  using  the  transfer 
matrix  method.  We  will  first  consider  the  behaviors  of  one¬ 
dimensional  conductance  at  normal  incidence,  and  suppress 
the  labels  and  cr. 

For  any  one-dimensional  channel  consisting  of  homoge¬ 
neous  dielectric  material,  @((u)  =  l  for  all  frequencies.  Equa¬ 
tion  (2)  then  leads  to  the  universal  quantized  thermal  con¬ 
ductance  Go(7’)  =  (6^).'’^’"^  For  any  one-dimensional 

channel  in  general,  we  define  its  normalized  thermal  conduc¬ 
tance  as  G{T)IGq{T). 

We  now  briefly  review  the  property  of  thermal  conduc¬ 
tance  suppression  of  a  single  photonic  crystal  at  normal 
incidence. The  suppression  is  particularly  effective  at  the 
ergodic  limit  of  kpja! {hc)>\,  where  the  thermal  conduc¬ 
tance  has  contribution  from  a  large  number  of  photonic 
bands.  In  this  limit,  the  statistical  properties  of  the  photonic 


band  structures  as  a  whole,  rather  than  the  detailed  properties 
of  individual  bands,  become  important.  From  Eq.  (2),  we 
thus  obtain 

7  f  =  r),  for  A  =  1.  (3) 

r-^OoCjQ^i  j  A^OoAJq 


Here  ?;<  1  denotes  the  proportion  of  frequencies  in  the  pho¬ 
tonic  bands.  Using  ergodic  theory,  we  showed  in  Ref.  10  that 
except  for  a  zero-measure  set  of  parameter  choices,  all  crys¬ 
tals  satisfy 


1  f" 


d(j)\  Re 


cos 


_i  -  2  H-  (^H-  l)cos  (f) 
^-1  . 


Re 


cos 


,,2-h  (g-H  l)cos  (j) 
^-1 


(4) 


where  g=0.5(nH- 1/n).  We  emphasize  that  77,  and  hence  the 
ergodic  limit  of  the  normalized  conductance,  is  independent 
of  the  layer  thicknesses  in  the  unit  cell. 

In  a  heterostructure  with  N  crystals  [Eig.  1(a)],  photonic 
band  gaps  are  located  at  different  frequencies  for  different 
crystals  [Eig.  1(b)].  To  determine  its  thermal  conductance  in 
the  ergodic  limit,  we  first  assume  that  thermal  transport  oc¬ 
curs  only  at  frequencies  that  fall  within  the  bands  of  all  crys¬ 
tals,  and  set  0(w)  =  O  for  all  other  frequencies.  We  then  as¬ 
sume  that  the  frequency  distributions  of  the  photonic  bands 
between  different  crystals  are  independent.  Erom  these  two 
assumptions,  the  total  proportion  of  frequencies  that  contrib¬ 
utes  to  thermal  conductance  is  then  rf,  with  77  given  by  Eq. 
(4).  Einally,  for  each  frequency  w  that  does  contribute,  the 
power  transmission  coefficient  0(w)  is  typically  less  than 
unity,  due  to  impedance  mismatch  between  different  crystals. 
Combining  all  these  considerations,  we  have 


G(T)  1 

lim — 7T  =  1™  T  I  du)&{u))<rf,  for  A  &  2.  (5) 

T^cdGq{T)  A^<»Adg 


Equation  (5)  is  the  main  result  of  this  paper:  the  thermal 
conductance  of  a  heterostructure  in  the  ergodic  limit  has  an 
upper  bound  rf  that  decreases  exponentially  with  the  num¬ 
ber  of  crystals  N.  Moreover,  such  an  upper  bound  is  indepen¬ 
dent  of  the  detailed  geometry  such  as  layer  thickness  and 
lattice  constant  but  instead  depends  only  on  the  refractive 
indices  of  the  layers. 

We  now  support  the  theoretical  analysis  above  with  ex¬ 
tensive  numerical  simulations.  We  first  verify  that  the  er¬ 
godic  limit  indeed  exists.  In  Eig.  2(a),  we  plot  the  tempera¬ 
ture  dependence  of  G{T)I  Gq{T)  for  various  heterostructures. 
G{T)IGoiT)  all  decrease  rapidly  at  small  T  and  converge  to 
specific  values  when  T> hcl  (k^a).  These  specific  values: 
lim7-^oo[G(7’)/Go(7’)],  which  we  refer  as  the  “high-T  conduc¬ 
tance,”  define  the  ergodic  limit  for  the  conductance  in  a 
given  structure.  Below,  we  will  only  consider  conductance  at 
this  ergodic  limit. 

In  a  heterostructure,  except  for  the  two  crystals  at  the 
ends,  all  intermediate  crystals  have  a  finite  number  of  peri¬ 
ods  Np.  At  frequencies  in  the  band  gaps  of  any  intermediate 
crystal,  photons  can  still  transmit  through  the  crystal,  result¬ 
ing  in  a  nonzero  0(&7).  Such  a  finite-size  effect,  however. 
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should  be  negligible  if  Np  is  sufficiently  large.  In  Fig.  2(b), 
we  consider  heterostructures  with  different  number  of  crys¬ 
tals  N.  For  each  N,  we  vary  Np.  The  high-T  conductance  all 
converges  when  ApS  10.  Thus,  with  a  choice  of  Np=  10,  all 
intermediate  crystals  can  essentially  be  treated  as  inhnite, 
and  thus  @((u)  =  0  unless  at  the  frequency  co  there  are  photon 
states  in  the  band  structures  for  all  crystals,  validating  the 
hrst  assumption  made  in  our  analysis.  Regarding  the  second 
assumption,  one  can  in  fact  prove  that  frequency  distribu¬ 
tions  of  photon  states  in  different  crystals  are  independent, 
except  for  a  zero-measure  set  of  parameter  choices,  using  an 
analysis  similar  to  Ref.  10. 

Equation  (5)  predicts  that  an  upper  bound  that  is  universal 
in  the  sense  that  it  is  independent  of  layer  thicknesses.  Nu¬ 
merical  calculations  establish  an  even  stronger  result.  In  Fig. 
2(c),  we  consider  a  total  of  15  000  structures,  with  the  num¬ 
ber  of  crystals  N  between  2  and  4.  In  these  structures,  the 
thickness  of  the  dielectric  layers  is  randomly  generated.  For 
the  vast  majority  of  structures,  the  high-T  conductance  them¬ 
selves  have  universal  values  that  depend  only  on  N  and  are 
independent  of  the  layer  thicknesses.  This  numerical  result 
provides  a  strong  support  of  our  statistical  analysis  that  leads 
to  Eq.  (5). 

Examining  Eig.  2(c),  we  also  see  few  cases  for  which  the 
high-T  conductance  of  a  structure  deviates  from  the  universal 
values.  [One  of  these  cases  is  also  shown  in  Eig.  2(a).]  Such 
deviations  occur  if  one  of  the  crystals  belongs  to  the  com¬ 
mensurate  cases  either  as  discussed  in  Ref.  10,  or  when  pho¬ 
tonic  bands  of  two  or  more  crystals  become  correlated  with 
each  other.  These  cases,  in  the  ergodic  limit,  form  the  zero- 
measure  set  as  discussed  above.  We  will  ignore  these  com- 


EIG.  3.  One-dimensional  high-T  conductance  for  photonic  crys¬ 
tal  heterostmctures  as  a  function  of  number  of  crystals  (N).  Data 
points  exist  only  at  integer  N.  N=0  corresponds  to  vacuum.  N=l 
corresponds  to  an  infinite  crystal.  N=2  corresponds  to  two  semi¬ 
infinite  crystals  placed  together.  For  Aa  3,  each  intermediate  crystal 
has  Ap=  10  number  of  periods.  The  solid  line  is  the  actual  conduc¬ 
tance.  The  dashed  line  is  the  corresponding  theoretical  upper  bound 
from  Eq.  (5),  where  q'=0.62. 

mensurate  cases  hereafter  since  their  occurrences  require 
control  of  layer  thicknesses  with  high  accuracy. 

We  now  provide  a  direct  numerical  confirmation  of  Eq. 
(5).  In  Eig.  3,  we  plot  the  high-T  conductance  of  heterostruc¬ 
tures  with  different  number  of  crystals  N.  We  consider  only 
the  incommensurate  cases  where  the  high-T  conductance 
takes  the  universal  values.  The  numerically  determined 
high-T  conductance  clearly  falls  below  the  exponentially  de¬ 
creasing  bound  as  set  by  Eq.  (5). 

Eor  three-dimensional  thermal  conductance,  using  Eqs. 
(1)  and  (2),  we  have'*’ 


kp 


0(CU,M,(7) 


(6) 


where  M  =  k||C/w.  At  the  ergodic  limit,  for  each  channel  la¬ 
beled  by  M,  we  repeat  the  same  statistical  arguments  above  to 
obtain 


G.uiT) 
hm - — 

7-^00  Gvac(T) 


fl  j  rA 

2  I  udu  ■  lim  Y  I  da)&((o,u,aj 

T=.'i,pJo 


—  2  I  [7j{u,cr)]^udu. 

(T=S,p  0 


(7) 


Here  rj{u,cr)  is  obtained  by  replacing  ^  in  Eq.  (4)  with 
^{u,a-)  =  0.5[nJ{ni,P)  +  {ni,P)/n^],  where  n^=\n^-u^,  tij, 
=  \T-m^,  and  P=\  and  n^  for  the  s  and  p  polarizations,  re¬ 
spectively. 

In  Eq.  (7),  the  equality  holds  only  when  N=\.  Also,  in 
deriving  Eq.  (7),  we  use  the  fact  that  in  the  ergodic  limit,  the 
states  with  m  >  1 ,  where  photons  are  evanescent  in  the 
vacuum  layers,  do  not  contribute.'"  Hence  the  range  of  inte¬ 
gration  on  u  is  between  0  and  1.  Einally,  Eq.  (7)  indicates  a 
universal  upper  bound  that  depends  only  on  the  index  and 
the  number  of  crystals.  Similar  to  the  case  of  normal  inci¬ 


dence,  numerical  results  (not  shown  here)  further  indicate 
that  3D  conductance  themselves  in  the  ergodic  limit  are  in 
fact  universal. 

We  now  confirm  Eq.  (7)  with  numerical  results.  Eor  the  s 
polarization,  the  band  gaps  persist  in  the  entire  range  of  0 
<M<1."''"  Hence  97(1/,  £7=  i)  <  1.  Erom  Eq.  (7),  the  upper 
bound  on  the  high-T  conductance  for  the  s  polarization  there¬ 
fore  decreases  exponentially  with  respect  to  the  number  of 
crystals  N,  as  supported  by  the  numerical  results  in  Eig.  4. 
Eor  the  p  polarization,  Up=\n^l (n^+l)  corresponds  to  the 
Brewster  angle  where  the  reflection  at  the  vacuum-dielectric 
interface  vanishes.'*’'®  At  u  =  up  there  is  no  photonic  band 
gap  at  any  frequency  and  Tj{up , a=  p)=l.  As  a  result,  the 
high-T  conductance  from  the  p  polarization,  and  hence  the 
total  3D  high-T  conductance,  show  more  gradual  decay  with 
respect  to  N  (Eig.  4). 

In  our  structure,  exponential  reduction  in  3D  thermal  con¬ 
ductance  can  be  accomplished  using  a  simple  angle  filter  that 
blocks  out  thermal  radiation  with  a  large  angle  of  incidence. 
As  a  simple  illustration,  restricting  the  integration  range  of  u 
in  Eq.  (6)  to  0  ^  m  s  0.5  is  sufficient  to  generate  the  exponen¬ 
tial  reduction  (Eig.  4).  Alternatively,  for  heterostmctures  be- 


113105-3 


BRIEF  REPORTS 


PHYSICAL  REVIEW  B  82,  1 13 105  (2010) 


FIG.  4.  (Color  online)  3D  high-T  thermal  conductance  versus 
number  of  crystals  (N)  in  the  structures.  Data  points  exist  only  at 
integer  N.  The  thick  black  lines  are  the  total  conductance.  The  thick 
green  (gray)  lines  are  the  total  conductance  with  an  angle  filter.  The 
thin  dark  blue  (dark  gray)  and  the  thin  red  (light  gray)  lines  are 
contributions  from  the  ^  and  p  polarizations,  respectively.  The  solid 
lines  are  direct  numerical  calculations  of  the  thermal  conductance. 
The  dashed  lines  are  the  theoretical  upper  bounds  from  Eq.  (7).  For 
Vs  3,  each  intermediate  crystal  has  Vp=10  number  of  periods. 

tween  two  semi-infinite  vacuum  regions,  we  can  achieve  ex¬ 
ponential  reduction  using  omnidirectional  reflectors  as  the 
intermediate  crystals. Finally,  Ref.  20  shows  that  the  Brew¬ 
ster  angle  can  be  made  imaginary  using  films  with  large  bi¬ 
refringence.  Multilayer  structures  consisting  of  such  birefrin- 
gent  films  and  vacuum  therefore  could  also  be  useful  for 
demonstrating  exponential  reduction  in  3D  thermal  conduc¬ 
tance. 

We  now  briefly  comment  on  the  experimental  feasibility 
in  demonstrating  our  predictions.  Multilayer  structures  can 
be  fabricated  by  a  variety  of  techniques.'*’^®  Our  theory  calls 
for  a  measurement  of  the  photonic  thermal  conductance  in 
these  structures.  For  this  purpose  we  notice  the  exciting  re¬ 
cent  experimental  developments  in  measuring  photonic  ther¬ 
mal  conductance. Here,  one  can  consider  a  set  up  similar 


to  Refs.  21  and  22,  where  thermal  contacts  are  made  to  the 
ends  of  the  structure,  and  thermal  conductance  is  measured 
by  recording  the  amount  of  power  transferred.  Our  results 
may  also  be  practically  significant  for  ultimate  thermal  insu¬ 
lation.  Currently,  an  effective  mechanism  for  thermal  insula¬ 
tion  utilizes  multiple  layers  of  metal  and  vacuum.^*’^"'  In  such 
a  structure,  thermal  transfer  is  incoherent.  Its  thermal  con¬ 
ductance  scales  as  1  /  where  is  the  number  of  metal 
layers. Our  results  here  indicate  a  qualitatively  different 
mechanism  for  achieving  thermal  insulation  using  dielectric 
systems,  which  may  present  new  opportunities.  For  example, 
the  structure  in  Fig.  1  allows  high  transmission  through  nar¬ 
row  band  of  frequencies,  which  may  be  important  for  com¬ 
munication  through  such  a  thermally  insulating  medium. 

We  end  by  highlighting  closely  related  works.  The  con¬ 
cept  of  using  a  single  phononic  crystal  structure  to  suppress 
coherent  phononic  thermal  conductance  was  discussed  in 
Ref.  25,  and  was  extended  to  the  photonic  case  in  Refs.  9  and 
10.  In  these  structures,  the  thermal  conductance  is  indepen¬ 
dent  of  the  channel  length.  Meanwhile,  mismatch  between 
two  different  phononic  band  structures  had  been  used  to 
achieve  very  small  interfacial  thermal  conductance. Our 
result,  showing  exponential  suppression  of  thermal  conduc¬ 
tance  using  a  large  number  of  crystals,  is  a  step  further  and 
points  to  a  new  regime  of  coherent  thermal  transport.  While 
we  have  considered  photonic  thermal  transport,  the  general 
concept  in  this  paper  may  be  applicable  to  other  thermal 
carriers  as  well  provided  that  the  coherent  effect  is  signifi¬ 
cant.  Finally,  in  analogy  to  electron  transport,  where  an  ex¬ 
ponential  reduction  in  electronic  conductance  with  respect  to 
distance  is  a  direct  signature  of  electron  localization,  our 
result  indicates  an  intriguing  notion  that  heat  can  be  localized 
as  well. 

This  work  was  supported  by  AFOSR-MURI  programs 
(Grants  No.  FA9550-08- 1-0407  and  No.  FA9550-09- 1-0704) 
and  by  AFRL.  W.T.L.  also  acknowledges  the  support  of  the 
Croucher  Foundation  in  Hong  Kong. 


*  wahtung  @  gmail.com 
'^jushen  @  gmail.com 
*shanhui  @  stanford.edu 

*  A.  Greiner  et  al,  Phys.  Rev.  Lett.  78,  1114  (1997). 

^L.  G.  C.  Rego  and  G.  Kirczenow,  Phys.  Rev.  Lett.  81,  232 
(1998). 

*M.  V.  Simkin  and  G.  D.  Mahan,  Phys.  Rev.  Lett.  84,  927  (2000). 
^K.  Schwab  et  al..  Nature  (London)  404,  974  (2000). 

^D.  G.  Cahill  et  al,  J.  Appl.  Phys.  93,  793  (2003). 

®T.  Yamamoto  et  al,  Phys.  Rev.  Lett.  92,  075502  (2004). 

’D.  R.  Schmidt  et  al.,  Phys.  Rev.  Lett.  93,  045901  (2004). 

*M.  Meschke  et  al.  Nature  (London)  444,  187  (2006). 

‘'W.  T.  Lau  et  al.,  Appl.  Phys.  Lett.  92,  103106  (2008). 

'“W.  T.  Lau  et  al.,  Phys.  Rev.  B  80,  155135  (2009). 

"W.  Kim  et  al..  Nanotoday  2,  40  (2007). 

*^C.  Kittel  and  H.  Kroemer,  Thermal  Physics  (W.  H.  Freeman,  San 
Francisco,  1980). 

'*X.  Wang  et  al,  Appl.  Phys.  Lett.  80,  4291  (2002). 

'^B.-S.  Song  etal.  Science  300,  1537  (2003). 


'^E.  Istrate  and  E.  H.  Sargent,  Rev.  Mod.  Phys.  78,  455  (2006). 

'^E.  D.  Palik,  Handbook  of  Optical  Constants  of  Solids  (Aca¬ 
demic,  New  York,  1985),  p.  554. 

'^A.  Yariv  and  P.  Yeh,  Photonics:  Optical  Electronics  in  Modern 
Communications,  6th  ed.  (Oxford  University  Press,  Oxford,  En¬ 
gland,  2007). 

'*Y.  Fink  et  al.  Science  282,  1679  (1998). 

'®J.  D.  Jackson,  Classical  Electrodynamics  (Wiley,  New  York, 
1998). 

20m.  F.  Weber  et  al.  Science  287,  2451  (2000). 

2'S.  Shen  et  al.  Nano  Lett.  9,  2909  (2009). 

22e.  Rousseau  et  al,  Nat.  Photonics  3,  514  (2009). 

2^0.  R.  Cunnington  and  C.  L.  Tien,  Prog.  Astronaut.  Aeronaut.  23, 
111  (1970). 

2"'C.  L.  Tien  and  G.  R.  Cunnington,  Adv.  Heat  Transfer  9,  349 
(1973). 

25  G.  Chen,  J.  Heat  Transfer  121,  945  (1999). 

2®V.  Narayanamurti  et  al,  Phys.  Rev.  Lett.  43,  2012  (1979). 

22h.  K.  Lyeo  and  D.  G.  Cahill,  Phys.  Rev.  B  73,  144301  (2006). 


113105-4 


